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ABSTRACT 

This paper describes a two-dimensional model developed to study the 
influence of convective flow on the concentration field ahead of a growing lamellar 
eutectic, when one phase projects out into the melt creating a stepped interface. The 
two dimensional convective flow field, which is periodic in the horizontal direction, 
was computed numerically using the software FLUENT. The velocity field generated 
due to the flow of melt over the steps was then incorporated into a finite difference 
scheme employed to solve the concentration field. The average interfacial composition 
was calculated and converted to lamellar spacing using the Jackson and Hunt mini- 
mum supercooling criterion. It was found that, a stepped interface is more sensitive 
to convection than a planar interface. 


1. INTRODUCTION 


The primary motivation for this work arises from experiments on eutectic 
solidification in the reduced gravity environment of space. Some space-grown samples 
exhibited marked structural differences from identically processed earth-grown sam- 
ples [1,2]. For example, when the MnBi/Bi eutectic was solidified in space the average 
fiber spacing A was half of what was obtained on earth [3]. The difference has been 
attributed to the absence of convective currents at low g as compared to conditions 
on earth. The application of a magnetic field to suppress convection showed the same 
effect on the MnBi/Bi microstructure as solidification in space [10]. 

As a result of the above experimental observations, theoretical models have been 
developed to determine the effect of convection on the spacing of eutectic alloys. 
Baskaran and Wilcox [4] and Chandrasekhar et al. [5] developed a two-dimensional 
model to study the effect of convection on the microstructures of lamellar eutectics. 
Caram and Wilcox [7] developed a three dimensional numerical model for the influ- 
ence of convection on rod-like eutectic microstructure. 

Decantation experiments during the solidification of MnBi/Bi eutectic estab- 
lished that the MnBi fibers project ahead of the Bi matrix [6]. In all the prior numer- 
ical models developed to predict the effect of convection on eutectic microstructure 
it had been assumed that the interface is planar. In the present work a lamellar 
stepped structure was studied. This was selected rather than a rod structure for 
computational ease; the diffusional domain of a lamellar eutectic is two dimensional 
rather than three dimensional as in the case of a rod eutectic. The Jackson and Hunt 
treatment [8] was assumed to be valid, and we performed calculations similar to those 
of ref. [5]. The resulting conclusions are expected to hold for a rod structure in a 
semi-quantitative sense. Details are given elsewhere [9]. 
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2. EQUATIONS 


As the adjacent phases of a binary eutectic grow, atoms are rejected 
into the melt due to partitioning. This partitioning creates a lateral composition gra- 
dient, which generates lateral diffusion in order to redistribute the two components 
across the interface and to sustain growth. As a result, the eutectic phase transfor- 
mation is influenced by mass transfer processes. In addition to diffusion, convective 
mass transfer generated by the flow along the solid liquid interface must be taken into 
account. 

The general steady state mass transfer equation for a binary mixture can be 
written as: 

DV 2 C - U ■ VC = 0 (1) 

where D is the diffusion coefficient, U is the velocity and C is the concentration of one 
of the components. To account for the influence of convection in eutectic growth, a 
well developed laminar shear flow was introduced in front of the solid-liquid interface, 
normal to both the growth direction of the lamellae and to the lamellae themselves. 
The velocity field generated by the flow of melt over the steps was computed nu- 
merically. The melt was assumed to possess constant material properties, including 
diffusivity D and kinematic viscosity v. The two dimensional flow field was assumed 
to be periodic in the horizontal direction. The governing fluid flow equations were 
written in a reference frame that moves vertically upward at the growth velocity. A 
no-slip condition was used for the tangential component of the velocity at the inter- 
face. The boundary at the top of the domain was placed at such a distance from the 
interface that it did not influence either the flow field or the concentration field near 
the interface. 

The fluid velocity components U„ in the y direction and U* in the x direction were 
calculated using a software package called FLUENT. FLUENT is a computer code 
created by Creare Incorporated for simulating a wide range of fluid flow problems. 
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It uses a finite difference numerical procedure to solve the Navier-Stokes equation. 
The domain to be investigated is divided into a finite number of cells and the partial 
differential equation is discretized over these units to produce a sequence of algebraic 
relations. An iterative scheme is then employed to solve for the equations and yield 
the velocity components U x and U y at each point in the finite difference grid over the 
problem domain. 

The velocities obtained as above were incorporated into the governing mass trans- 
fer equation, which was solved using finite difference techniques. The governing partial 
differential equation for mass transfer was simplified as: 

d 2 C d 2 C Uy dC U x 3C _ 0 / 9) 

dx 2 + dy 2 D dy D dx 

The boundary conditions used for determining the solute concentration in the 
melt ahead of the solid liquid interface are as follows. 

The composition remains at the eutectic far from the interface: 

C = C e at y = oo (3) 


Due to conservation of matter at the interface, the flux of solute rejected at the 
freezing interface is equal to the solute diffusing into the melt away from the interface. 
This condition can be mathematically represented in the following form: 

For a phase: 

D (f )v=0 = -V (C, - Cf) (4) 

For /? phase, with a stepped interface: 

D = -V (C, - Cf) (5) 

There is no side growth, i.e. along the side of the step there is no diffusion in the x 
direction: 


( dx ^° to 


= 0 
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( 6 ) 



The concentration field is periodic owing to the periodicity of the lamellar structure: 


C , (x=0,y) — C , (x=A,y) — C(x=n\,\i) n — 1>2, •• (7) 

The far field boundary condition, eq. 3, was tried for a stepped interface for 
increasing values of y until the solution was independent of the domain size. It was 
found in the initial computations that the computed average interface concentration 
with a far field condition at y=3A/4 differed by less than 1% from that computed 
with y=A/2. Hence in the subsequent computations we used C=C e at y= 3A/4. 

The solid solubility was assumed to be negligible, so C* =0 and Cjj=l. The de- 
gree of non-planarity of the interface was altered by giving various values to «, which 
is the height by which the alternate lamellae project into the melt. If we define k — 
cA, then c was chosen in the range 0 to 1, because Chandrasekhar [6] reported that 
the length-to-diameter ratio of the protruding MnBi rods averaged one. 

In order to determine the effect of convection on the lamellar spacing we min- 
imized the average interfacial undercooling AT. The treatment is similar to that in 
Chandrasekhar and Wilcox [5]. To facilitate the numerical computations we non- 
dimensionalized equation (2) to the following form, with the symbols defined in the 
table of nomenclature: 


d 2 C d 2 C . 

h 1- A 

dX 2 dY 2 


— -ry = o 
dY dx 


( 8 ) 


Since compositional changes occur over distances on the order of A } the dimensionless 
distances were defined by X=x/A and Y=y/A. Note that A= XV/ D is a lamellar 


spacing Peclet number and T = G U A 2 /D is a convective Peclet number. Letting 
S a + 2 Sfi = S , the non-dimensional forms of the boundary conditions are: 

at Y = 3/4 

C = C t 

(9) 

at Y = 0 

0 < X < Sa/ A, S/X < X < 1, 

W=- ACi (10) 

0 < Y < «, 

dC 

X = SJ A or S/X, 0 

(11) 
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SJX < X < 6/ A, 




( 12 ) 


For typical growth conditions on earth A is on the order of 0.1 and T is of order 100. 

3. Results and discussion 

Figure 3 shows a typical result, the the iso-concentration contours for a eutec- 
tic of C e =0.3 with a step height of A/8 growing without the disturbing influence of 
convection. Convection distorts the contours as shown in figure 4 for a convective 
Peclet number T = 40. The computed concentration fields in the melt were converted 
to lamellar spacing by calculating the average interfacial composition Ci and mini- 
mizing the total interfacial undercooling AT. Following Chandrasekhar et al. [5] the 
deviation parameters are defined as follows: 

A. = ( Ci. -C,), - ( Ci. -a) (!3) 

A„ = ( Ci, -c«)„ - ( Ci, -C') (H) 

where A a is the difference between the deviation from the eutectic composition with- 
out convection (G u = 0) and that with convection (G u > 0) averaged over the a 
phase. Here A p is the difference between the deviation from the eutectic composition 
without convection (G u = 0) and that with convection (G u > 0) averaged over the j3 
phase. The subscript o indicates the values without convection. 

The deviation from eutectic composition without convection was found to be 
proportional to A = A V/D: 

( Ci a -C e )o = A a A (15) 

( Ci, ~C e ) 0 = A 0 A (16) 

The analytical solution of Jackson and Hunt [8] gives the values of A a and A 0 for a 
planar interface and no convection. From Jackson and Hunt [8] A = 0.121 for C e =0.3 
and A = 0.135 for C e =0.5. From our numerical results, A = 0.119 for C e =0.3 and A 
= 0.136 for C e = 0.5 for the planar interface. For a stepped interface with C e =0.3 
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and #c=A/5 we found A = 0.129. For C e =0.3 and /c= A/4, A = 0.134. 

Following Baskaran and Wilcox [4] the spacing between the lamellae is given by: 



Differentiating f with respect to T we find df/dr and substitute this into eq. (17). 
Polynomial fits for A/A 0 as a function of T 0 are shown in figures 5 to 7. 

For C e =0.1 and a planar interface: 

A/A 0 = 1 + 0.0016 r 0 + 1.61 X 10- 5 T 2 C 
For C e = 0.1 and «= A/8: 

A/A 0 = 1 + 0.0019 r 0 + 2.28 X 10" 5 
For C e = 0.5 and a planar interface: 


A/A 0 = 1 + 0.0073 T 0 - 0.00015 T 2 0 
For C e = 0.5 and /c=A/8: 

A/A 0 = 1 + 0.0114 r o - 0.0003 
For C e = 0.3 and a planar interface: 


A/A e = 1 + 0.00046 r o + 9.80 X 10“ 5 T 2 0 
For C e = 0.3 and /c=A/5: 


x/x 0 = i + o.oo76 r e - o.ooo5 r* 


4. Conclusions 


Based upon the results obtained from the numerical analysis the fol- 
lowing conclusions can be made: 
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1 . Increasing the intensity of convection increases the spacing between the lamellae. 

2. For the same convective Peclet number, the increase in lamellar spacing of a 
eutectic growing with a stepped interface is much higher than the increase in 
lamellar spacing observed for a eutectic growing with a planar interface. That 
is, the microstructure of stepped interface is more sensitive to convection than 
a planar interface. 

3. The increase in lamellar spacing observed for the same value of convective P eclet 
number is more for the eutectic growing with a higher k, i.e. more non-planar. 
That is, with increasing irregularity of the interface the eutectic experiences an 
increased change in the lamellar spacing due to convection. 

4. The effect of convection is negligible if the flow velocity is very small in magni- 
tude, i.e. T is small (<1). 
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Nomenclature 

A a = constant in equation 15 

Ap = constant in equation 16 

C e = Eutectic mass fraction 

Cj = Interfacial melt composition 

Ci = Average interfacial melt composition 

C“ = Composition of the alpha phase (assumed 0 here) 

Cf = Composition of the beta phase (assumed 1 here) 
f = A/A 
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G u = Gradient of convective velocity in the y direction [s -1 ] 
S a = Half width of a phase [m] 

T = Temperature [K] 

AT = Total interfacial undercooling [K] 

U = Convective velocity [m/s] 

V = Growth rate [m/s] 

x = Distance along interface from beginning of a phase [m] 

X = Dimensionless distance along interface (x/A) 
y = Distance into the melt from interface [m] 

Y = Dimensionless distance into the melt from interface (y/ A) 
a = The A-rich phase of the binary eutectic 

/? = The B-rich phase of the binary eutectic 
S = S a *t* 2S/3 

A = ( Ci a ~C e ) 0 - ( Ci a — C e ), the deviation parameter 
r = G U A 2 /D, Peclet number for convection 
r 0 = G U A 2 /D 

k = The height of the step, as a fraction of A 
A = Lamellar spacing [m] 

A = A V /D Lamellar spacing based Peclet number 
H = Viscosity of the melt ( 4 X 10 -3 kg/m.s used here) 
p = Density of the melt ( 800 kg/ m 3 used here) 
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